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ABSTRACT 

Using a modification of the Zel'dovich approximation adapted to hot dark 
matter, the accretion of such matter onto moving cosmic string loops is studied. 
It is shown that a large number of 1O 12 M0 nonlinear objects will be produced by 
a redshift of z = 4. These objects could be the hosts of high redshift quasars. 
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1. Introduction 



Observation of high redshift objects has emerged as a powerful tool for testing 
theories of structure formation. For example, in inflationary Universe models, the 
epoch when the first massive nonlinear structures (which could be the hosts of 
quasars or primordial galaxies) form is a sensitive function of the fraction of hot 
dark matter. Recent data on the abundance of damped Lyman alpha absorption 
systems (DLAS) 1 ) and on the quasar abundance 2 ) has provided tight constraints on 
inflation-based models of structure formation with adiabatic density fluctuations 
and mixed dark matter. 

Searches for high redshift quasars have been going on for some time 2 ). The 
quasar luminosity function is observed to rise sharply as a function of redshift 
z until z ~ 2.5. According to recent results from the Palomar grism survey by 
Schmidt et al. (1995) 2 ), it peaks in the redshift interval z e [1.7, 2.7] and declines at 
higher redshifts. Irwin et al. 2 ) on the other hand find that the luminosity function 
is constant up to redshifts of about 4. Quasars (QSO) are extremely luminous, 
and it is generally assumed that they are powered by accretion onto black holes. 
It is possible to estimate the mass of the host galaxy of the quasar as a function of 
its luminosity, assuming that the quasar luminosity corresponds to the Eddington 
luminosity of the black hole. For a quasar of absolute blue magnitude Mb = —26, 
the host galaxy mass can be estimated as 3 ) 

M G = cilO 12 M , (1.1) 

where c\ is a constant which contains the uncertainties in relating blue magnitude 
to bolometric magnitude of quasars, in the baryon fraction of the Universe and 
in the fraction of baryons in the host galaxy able to form the compact central 
object (taken to be 10~ 2 ). The best estimate for c\ is about 1. Models of structure 
formation have to pass the test of producing enough early objects of sufficiently 
large mass to host the observed quasars. 
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Damped Lyman alpha systems provide potentially even more powerful con- 
straints on structure formation models. Evidence is mounting 4 ) that the absorption 
line systems observed in the spectrum of distant quasars are due to progenitors of 
typical galaxies. Based on the number density of absorption lines per frequency 
interval and on the column density calculated from individual absorption lines, 
the fraction of fl in bound neutral gas (denoted by fl g ) can be estimated. Recent 
observational results 1 ) indicate that 

Q g (z) > 1(T 3 (1.2) 

for ze [1,3], with the highest value in fact taken on at z = 3 ! Note that the above 
corresponds to a value Q in bound matter which is larger by a factor of f^ 1 , where 
is the fraction of bound matter which is baryonic. In a Q = 1 cosmology, the 
value of ft, is expected to be of the order 10 _1 . 

The constraints coming from observation of high redshift objects for inflation- 
based models of structure formation were studied by several groups. In Ref. 7, it 
was shown that the high redshift quasar abundance is compatible with an unbiased 
cold dark matter (CDM) model , but that the theory predicts an exponential 
decrease of the quasar abundance for z > 5. Recently 8 ), this analysis has been 
extended to mixed dark matter (MDM) models. The abundance of damped Lyman 
alpha systems was used in Refs. 3, 5 & 6 to further constrain MDM models. It was 
found that models with a fraction of hot dark matter (HDM) exceeding Vt v = 0.2 
do not predict enough nonlinear structures at high redshift in order to be able to 
explain the data. 

However, there exists a viable class of alternative theories of structure forma- 
tion, the cosmic string (CS) models. In these models, density fluctuations are gen- 
erated by topologically stable strings of trapped energy density (one-dimensional 
topological defects in a relativistic field theory describing matter), instead of origi- 
nating as the result of quantum fluctuations during an early period of exponential 

c) A model with flat spectrum (n = 1), ft = 1, and vanishing cosmological constant A. 
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expansion of the Universe. The main purpose of this paper is to present a prelim- 
inary study of the constraints on cosmic string models which can be derived using 
the QSO and DLAS abundances. The model we are most interested in is a cosmic 
string theory in a spatially flat (Q = 1) Universe with only HDM and baryons. 
This model is briefly reviewed in Section 2. 

Our main result is that for reasonable values of the parameters of the model, 
the CS and HDM theory is compatible with the present observational constraints 
on the quasar abundance of the Palomar grism survey by Schmidt et al. 2 ), which 
are plotted in Figure 1. We will compare these observations for the co moving space 
density of quasars brighter than Mb = —26 with the abundance of objects of mass 
greater than Mq (given in (1.1)) in the cosmic string and hot dark matter model.. 
We also comment on the implications of the DLAS abundance (as given by (1.2)) 
for cosmic strings. 

Our result is quite nontrivial. The reason why inflationary MDM models with a 
large fraction of HDM do not predict a sufficient abundance of nonlinear objects at 
high redshifts is that the spectrum of density perturbations is suppressed at small 
wavelengths (those which first become nonlinear in the standard CDM inflationary 
model) by neutrino free streaming 9 ). The primordial power spectrum in a cosmic 
string theory is scale invariant 10 \ as it is in inflation-based models. The reason 
why a CS and HDM model is viable at all is that cosmic strings survive free 
streaming 11 ^. Since the strings are long-lived seed perturbations (as opposed to 
adiabatic dark matter fluctuations), accretion of dark matter on small scales - 
wavelength A < Xf^, where \™ ax is the maximal neutrino free streaming scale 
(whose value is given later) - is delayed but not prevented. Thus, the spectrum 
of density perturbations is not exponentially cut off below \y ax as it is for an 
inflationary HDM theory. It is 12 -*, however, suppressed by a power of A/Aj" aa; 
compared to that of an inflationary CDM model. It also has a smaller amplitude 
than the MDM model with VL V = 0.2. Hence, it would seem that the CS and HDM 
model would be unable to explain the abundance of high redshift QSO and DLAS. 
However, as explained in the following paragraph, the above reasoning misses a 
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crucial point. 

In the cosmic string theory - in contrast to inflation-based models - the density 
field is not a Gaussian random field. There are localized high density peaks even 
when the average density contrast is small. Hence, knowledge of the density power 
spectrum is insufficient to calculate the number density of nonlinear objects. 

In particular, cosmic string loops seed large amplitude local density contrasts. 
In this paper, we study the accretion of hot dark matter onto moving string loops 
and use the results to compute the number density of high redshift objects as 
a function of a parameter v which determines the number density of loops in 
the scaling solution (see Section 2). We demonstrate that for realistic values of 
u, the number of massive nonlinear objects at redshifts < 4 satisfies the recent 
observational constraint of quasar abundances (see Figure 1). We also comment 
on the implications of (1.2). 

The next section of this paper contains a brief review of the CS and HDM 
theory of structure formation. In Section 3 we summarize the methods used: 
a Zel'dovich approximation 13 ) technique modified 14 ) to HDM and its adaptation 
to moving seed perturbations 15 ) . Section 4 contains the main calculations, and 
in Section 5 we discuss the results. Units in which h = c = ks = 1 are used 
throughout, and a Hubble constant of H = 50 /i5okms _1 Mpc _1 and a redshift at 
equal matter and radiation of z eq = 5750 £lh§Q are used. 

2. Brief Review of the Cosmic String and Hot Dark Matter Theory 

Cosmic strings 16 ) are one-dimensional topological defects which are predicted 
in many relativistic field theories describing matter. In such theories, a network of 
strings forms during a phase transition in the very early Universe. These strings 
are characterized by a mass per unit length /i (which in principle is the only free 
parameter in a cosmic string model) which determines their gravitational effects. 
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After the time of formation, the network of strings rapidly approaches a "scal- 
ing solution" , a distribution of defects whose statistical properties are time inde- 
pendent when lengths are rescaled by dividing by the Hubble radius. The mean 
separation of the strings increases by having strings interconnect and chop off loops. 
The scaling solution of the string network is characterized by a fixed number N of 
long strings crossing each Hubble volume at any time t, and the presence of loops 
with a distribution 

n(l,t) = vl~ 2 r 2 , (2.1) 

where / is the length of the loop and v is a constant. The quantity n(l,t)dt gives 
the number per unit physical volume at time t of loops with lengths in the interval 
between / and I + dl. The loops are remnants of string intercommutations at times 
t' < t. Loops oscillate and decay slowly by emitting gravitational radiation. Hence, 
there is a lower cutoff value of / for the distribution (2.1) given by 

Imin ~ Gflt , (2.2) 

G being Newton's constant. Below l m i n , the distribution n(l,t) becomes constant. 
We shall not discuss this point since it will not be relevant to our computations. 

In principle, the properties of the scaling solution and hence also the value 
of the constant v is calculable, albeit only numerically. In practice, however, the 
dynamics of the defect network is quite complicated and the numerical resolution 
inadequate to solve this problem. Thus, we must treat v as an undetermined 
constant. There are two more such constants which are denoted by a and (3. 
The first constant determines the mean radius Rf of a string loop at the time of 
formation 

R f (t) = at, (2.3) 
the second relates the mean radius R to the length / of a loop: 

l = (3R. (2.4) 
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Numerical simulations 17 ) indicate that a < 10~ 2 and (3 ~ 10. They also indicate 
that N ~ 10. From these values it follows that - unless v is extremely large - 
most of the mass of the string network resides in long strings (where long strings 
are defined operationally as strings which are not loops with radius smaller than 
the Hubble radius). 

In the cosmic string theory there are two basic mechanisms which seed struc- 
tures, loops and wakes. String loops act gravitationally almost like point mass 
objects when viewed from distances larger than R. Long straight strings, on the 
other hand, lead to planar over-dense regions called wakes 18 ). On distance scales 
smaller than its curvature radius, the local gravitational attraction of a string 
vanishes. However, space perpendicular to the string is conical with deficit angle 
87rG/i. Hence, a string moving with transverse velocity v s will impart a velocity 
perturbation 

Sv = AtiG^v s ^(v s ) (2.5) 

towards the plane behind the string, where 7(f) is the relativistic 7 factor. This 
develops into a planar over-dense region behind the string, the wake. 

Since most of the mass in the string network is in long strings, the wake mech- 
anism will be responsible for the formation of most of the present structure in the 
Universe. The thickest and most numerous wakes are those created at the time of 
equal matter and radiation t eq 19 \ The cosmic string theory hence predicts a dis- 
tinguished scale and topology of the large-scale structure in encouraging agreement 
with the data from the Cf A galaxy redshift survey 20 ) . 

The scaling for the defect network leads to a scale invariant spectrum of density 
perturbations, which in turn leads to a "scale- invariant" (n = 1) spectrum of 
microwave anisotropics 21 ). Normalizing the model by the COBE results gives 22 ) 

G/i~10~ 6 . (2.6) 
The accretion of hot dark matter onto cosmic string wakes was considered in 
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detail in Ref. 14. It was found that the first comoving scale to go nonlinear about 
a wake caused by a string at time t eq is 

qmaxiteq) = V eq Z eq t eq ~ V eq ■ 50h^ MpC (2.7) 

where v eq is the mean hot dark matter velocity at t eq . In a Q = 1 Universe, v eq 
is about 0.1. Hence, the distance q max is in good agreement with the observed 
thickness of the CfA galaxy sheets 20 -*. 

Note, however, that the first dark matter nonlinearities about wakes form only 
at late times, at a redshift 14 ) 

Zmax = ^-G{lV s '~i(v s )v~ q l Z eq , (2.8) 

which for v s l(v s ) ~ 1 and for the value of Gfj, from (2.6) corresponds to a redshift of 
about 1. Before this redshift, no nonlinearities form as a consequence of accretion 
onto a single uniform wake. 

Thus, in the cosmic string and hot dark matter theory, a different mecha- 
nism is required in order to explain the origin of high redshift objects. Possible 
mechanisms related to wakes are early structure formation at the crossing sites of 
different wakes 23 -* , small-scale structure of the strings giving rise to wakes 24 ' 25 ' 26 -* , 
and inhomogeneities inside of wakes 26 ). In this paper, however, we will explore a 
different mechanism, namely the accretion of hot dark matter onto loops. 

In earlier work, 11 ) the accretion of hot dark matter onto static cosmic string 
loops was studied. It was found that in spite of free streaming, the nonlinear struc- 
ture seeded by a point mass grows from inside out, and that the first nonlinearities 
form early on (accretion onto string filaments proceeds similarly 25 -*). In the con- 
text of the "old" cosmic string scenario (wakes unimportant), this mechanism was 
used in Ref. 27 to derive the mass function of galaxies. Since loop accretion leads 
to nonlinear structures at high redshift, we will now investigate this mechanism 
in detail to see whether enough high redshift massive objects to statisfy the QSO 
constraints and (1.2) form. 
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3. Modified Zel'dovich Approximation 

We will use the Zel'dovich approximation 13 ) and modifications thereof to study 
the accretion of hot dark matter onto moving string loops. The Zel'dovich approx- 
imation is a first order Lagrangian perturbation theory technique in which the 
time evolution of the comoving displacement ip of a dark matter particle from the 
location of the seed perturbation is studied. 

The physical distance of a dark matter particle from the center of the cosmic 
string loop is 

h(q,t)=a(t)[q-^(q,t)]. (3.1) 

The scale factor a(t) is normalized such that a(to) = 1, where to is the present time. 
The Zel'dovich approximation is based on combining the Newtonian equation for 
h 

h = -V h $ (3.2) 

with the Poisson equation for the Newtonian gravitational potential <£> and lin- 
earizing in ip. For a pointlike seed mass of magnitude m located at the comoving 
position q' — the resulting equation is 

^ + 2^ + 3^ = ^- (3.3) 

This equation describes how as a consequence of the seed mass, the motion 
of the dark matter particles away from the seed (driven by the expansion of the 
Universe) is gradually slowed down. If the seed perturbation is created at time t% 
and the dark matter is cold, then the appropriate initial conditions for ip are 

il>{q,U) =j){q,ti) = 0, (3.4) 

leading - for a(t) ~ t 2 / 3 appropriate in the matter dominated epoch t > t eq - to 
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the solution 



(3.5) 



As formulated above, the Zel'dovich approximation only works for cold dark 
matter, particles with negligible thermal velocities. The theory of interest to us, 
however, is based on hot dark matter. Luckily, the Zel'dovich approximation can 
be modified for HDM 14 ). HDM particles have large thermal velocities. At time t, 
the free streaming length in comoving coordinates is 

Xj(t) = v(t)z(t)t , (3.6) 

where v(t) ~ z(t) is the hot dark matter velocity. The length Xj(t) is the mean 
distance an HDM particle will move in one expansion time. Free streaming erases 
density perturbations on scales q < Xj(t), a scale which decreases as t -1 / 3 as t 
increases. 

A simple prescription ^ for taking into account free streaming in the Zel'dovich 
approximation is to - for a fixed comoving scale q - only let the perturbation start 
to develop at time t s (q) when 

Xj(t s (q)) = q, (3.7) 
i.e., replace the initial conditions (3.4) by 

^(q,t s (q))=^(q,t s (q)) = (3.8) 

with 

t s (q) = max{U,t B (q)} . (3.9) 

d) Note that this prescription has been shown 14 ) to give good agreement with an analysis 
obtained by tracking the full phase space distribution of HDM particles by means of the 
linearized collisionlcss Boltzmann equation. 
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The result for the comoving displacement ip{q) then becomes 

^=^ii)\mT- (310) 

We can now define the mass that has gone nonlinear about a seed perturbation 
as the rest mass inside of the shell which is "turning around", i.e. for which 

h(q,t) = 0. (3.11) 

This yields an equation 

q = 2^(q,t) (3.12) 

for the scale q n i(t) which is turning around at time t. For cold dark matter, Eq. 
(3.5) can be combined with Eq. (3.12) to obtain q n i(t) as well as the corresponding 
mass 

2 /^ 2/3 



M CDM (t) = -m ^-j . (3.13) 

Note the similarity of this result to what can be obtained in linear perturbation 
theory: (t/t^) 2 / 3 is precisely the growth factor of linear cosmo logical perturbations 
on small scales. 

For HDM, Eqs. (3.8), (3.10) and (3.12) can be combined to yield 

M »^ = w,mS < 3 - 14 > 

with 

M e , = ?%. (3.15) 

A further complication is due to the finite velocity of the cosmic string loops. 
This implies that we must extend the Zel'dovich approximation technique to mov- 
ing sources. For CDM, this was done by Bertschinger 15 ) with the interesting result 
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that there is to a first approximation no change in the total mass accreted. The 
suppression of the growth of perturbations in the direction perpendicular to the 
direction of motion of the seed mass is cancelled by the larger length of the non- 
linear region in the direction of motion. For HDM, however, there will be a net 
suppression of growth if the seed mass is moving. It will be important for us to 
take this effect into account. 

Given a moving point source, the basic Zel'dovich approximation equation (3.3) 
becomes a vector equation 



a ■ a Grain — q') 
i> + 2-4, + 3-^ = o, ~ }' (3.16) 
— a— a— cr|g — q'\° 



where q'(t) indicates the comoving position of the source. Without loss of generality 
we can take the source to move along the z axis with initial velocity V{ at time ti, 
so that 15 ) 

?'(*) = toft ^1 - (j^J ' j e z z(ti) , (3.17) 

e z being the unit vector along the z axis. In the matter dominated epoch Equation 
(3.16) can be solved exactly 15 ) (taking q y = without loss of generality) 



V>(<Z,t) = b(t)di[f x (q,t)e x + f z (q,t)e z ] (3.18) 



with 



b(t) 



di = 3viUz(ti) 

1 Gm a(t) (3.19) 



15vfti a(U) 



and where f x and f z are known functions of q, t and di which at late times become 
independent of time and contain the information about the geometry of the pattern 
of accretion onto the moving loop. In particular, the transverse displacement at 
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late times approaches 

'•<*> = it 



2 

Qx 



-^{Roo — Ri) + 



q x q z ~ 







Ridi_ 







(3.20) 



. X 2 
In the above, 

i2(9,0 = [?2 + (?*-^(0) 2 ] 1/2 , (3-21) 

i?i = R(q, ti), Roo = R(q, too), and k(q) is the factor by which the accretion at q is 
suppressed due to the motion of the source. 

It is easy to check that for q 2 3> df the factor k(q) tends to 1 (for q 2 3> q 2 ) and 
that the result for ip(q, t) from (3.18) goes to the result of (3.5). For rapidly moving 
loops we are interested in the opposite limit df 3> q 2 . In this case, evaluated at 
q z = 0, the suppression factor becomes 

k(q x )^2^-. (3.22) 

In this case, the "turn-around" condition 

q x = 2i/j x (3.23) 

for transverse accretion yields 



2 . . , , (to \ ft 



2/3 

For a cosmic string loop formed at time U whose mass is (see (2.3) and (2.4)) 

m = afifiti (3.25) 

this gives 

/fi \ 1/2 

q x ( ti ,t) = l^aP) (G^v-^toz-^it) . (3.26) 

From this result we can immediately recover Bertschinger's result 15 -* that accretion 
of CDM onto a moving loop is not suppressed: the total accreted mass M(t) is 
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proportional to 



M(t) ~ <£(t)dip c 



(3.27) 



where p c is the background comoving energy density. The factors of V{ evidently 
cancel! 

We are interested in the accretion of HDM onto moving loops. Provided that 
the turn-around distance q x of (3.26) exceeds the initial comoving free streaming 
distance, i.e. 



then the above analysis can also be applied to HDM. We will check this condition 
in our calculations in the following section. 



At this point, we are able to use the methods described in the previous section 
to compute the number density uq{> M\,t) in nonlinear objects heavier than Mi, 
and the fraction Q n i of the critical density in such objects at high redshift z(t) in 
the cosmic string and hot dark matter model. By considering only the accretion 
onto string loops we will be underestimating these quantities. 

It can be shown 27 ) that in an HDM model string loops accrete matter indepen- 
dently, at least before the large-scale structure turns nonlinear at the redshift given 
by (2.8) which is about 1, i.e., later than the times of interest in this paper. Hence, 
the number density n(l,t) of loops of length I given by (2.1) can be combined 
with the mass M{1) accreted by an individual loop (which follows from (3.13) and 
(3.14)) to give the mass function n(M,t). Here, n(M,t)dM is the number density 
of objects with mass in the interval between M and M + dM at time t. This in 



q x (U,t) > Xj(ti) 



(3.28) 



4. Computations 
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turn determines the comoving density in objects of mass M > Mi, 

M 2 



n G (> M h t) = z~ 3 (t) J dMn(M) (4.1) 



Mi 

and the fraction of the critical density in objects of mass greater than Mi, VL n [(t) : 

M 2 

n nl (t) = 6nGt 2 J dMMn(M) . (4.2) 
Mi 

If we want to compare with the observational results from QSO counts we must 
use appropriate integration limits Mi and M2 in (4.1) and (4.2). The mass Mi is 
the mass limit corresponding to the limiting QSO luminosity of the observational 
sample given in (1.1), Mi = ci10 12 Mq. For comparison with damped Lyman alpha 
systems, the lower cut-off mass Mi is somewhat smaller, but we will not need the 
exact value. For large loop radii, the approximation of treating the loop as a point 
mass breaks down. This will lead to a time dependent upper mass cutoff M2(t). 
A rough criterion for M2(t) can be obtained by demanding that the mass M(t) 
accreted onto a loop exceed the mass in a sphere of radius equal to the loop radius 

Using the CDM mass formula (3.13) with initial mass (3.25) yields the estimate 
M 2 (t)^|(0 1/2 (^) 3 / 2 |^)" 3/2 . (4-4) 



S-IO^M©, (4.5) 



Since 

to 
G 

the value of M2 is much greater than both Mi and the maximal neutrino Jeans 
mass M'(t), the largest mass which is affected by free streaming. This justifies the 
use of the CDM mass formula (3.13). 
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As indicated above, there is another mass which is crucial for the computation 
of ng(t) and Q n i(t), namely M'{t). The easiest way to determine M'(t) is to ask 
for the value of M for which the mass formulas (3.13) and (3.14) intersect: 



-m , 

5 \U 



2/3 



2\ 3 m 3 



5/ M% \t eq 



(4.6) 



Using the expression (3.25) for m and taking account of (3.15), one finds 



M'(t) 



2 /5 
5 V9 



1/4 



(a^) 3 / 4 4V*)- 3/ V /4 |. 



(4.7) 



from which it indeed follows that - at least for the redshifts we are interested in - 
M'(t) is smaller than M 2 (f). 

As a self-consistency check we note that M\ is larger than the mass accreted 
onto a loop created at t eq . Hence, it is consistent to restrict our attention to the 
matter epoch t > t eq . 

Since the functional form of M(m) and hence M(l) changes at M — M' , the 
functional forms of the mass function n(M) will be different above and below M' . 
Thus 

~ M'(t) M 2 (t) 

dMn H (M)+ J dMnc(M) . (4.8) 



n G (> M 1 ,t)=z~\t) 
Similarly for Q n i(t), 



Mi 



M'(t) 



Unlit) 



M'(t) 



M 2 (t) 



dMn H (M)M+ J dMn c (M)M 



Mi 



M'(t) 



(4.9) 



where njj{M) and nc{M) respectively refer to the mass functions computed with 
the HDM and CDM mass formulas (3.13) and (3.14) respectively. Combining (2.1) 



1G 



with (3.13) and (3.14) we obtain 



24 ifi 

n c (M,t) = v—(a[3) 2 ^ I (4.10) 

and 

*') = 4TP^jP' (4 ' u) 

As a consistency check, we note that for M = M'(t) the above two expressions 
coincide. 

By inspection of (4.8-4.11) it is clear that the integrals for no(t) and £l n l(t) 
are dominated at M — M'(t) and that a reasonable approximation to (4.8) will be 

n G (> Mi, t) ~ n c (M', ()lV 3 (i) ~ I^)2_^^-3 (t) ? (412 ) 
which gives 

n G {> Mi,t) ~ 5K«/3)- 1/4 (G^) 3 / 4 (^) 9 / 4 ^ 3 / 4 (t)t 3 

Ue 9 (4.13) 

~ 2 • 10- 4 ^- 3 / 4 (i)/i5/ 2 (/i5 1 Mpc)- 3 

when inserting t> e(? = 0.1, G// = 10~ 6 and the values of a = 10~ 2 and (3 — 10 from 
Section 2. This result is plotted in Figure 1. Similarly, (4.9) can be approximated 
by 

L 

which gives 



tord(t) ~ n c (M', t)M' 2 6nGt 2 ~ j 1 ^ 6rrGt 2 , (4.14) 



^"4(^ ,! « I,! (^) i(! *r J,! 
~ io _1 /ii i/^(*)" 3/2 

for the same values of the parameters as above. 



(4.15) 
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At this point it looks like the cosmic string and hot dark matter model will 
produce too many quasar host galaxies at redshifts z ~ 5. However, so far loop 
velocities have been neglected. 

Loop velocities can be taken into account by incorporating the condition (3.28). 
Loops which do not satisfy this criterion will not be able to accrete much mass. 
Note that as U increases, q x (U) remains constant whereas Xj(U) decreases. Hence, 
the condition (3.28) corresponds to a low mass cutoff M c in the integrals (4.1) and 
(4.2). Using (3.6) and (3.26) it follows that the inequality (3.28) becomes 

/fi \ 1//2 

Uz{ti)v{U) < (Jaf3) (Gpi^v-^tozit)- 1 / 2 , (4.16) 

where v{tj) is the hot dark matter velocity at time U and V{ is the loop velocity at 
the same time. After some algebra, (4.16) translates to 

z{ti) < ^aPGfizl^v-fzit)- 1 = z(t c ) (4.17) 

where v eq is the HDM velocity at t eq . The corresponding cutoff mass is 

M c = M(t c ) = laf3G^z{t c )- l l 2 z{t)- l % . (4.18) 
5 G 

It is easy to check that M c is larger than M', hence justifying the use of the CDM 
mass formula to obtain (4.18). This, however, also leads to a significant reduction 
in the values of no(t) and Q n i. 

The effect of loop velocities is thus to replace the estimate (4.12) by 
n G (> M u t) ~ n c (M c ,t)M c z- 3 (t) ~ 4z/(«/3) 1/2 (G/i) 3 /V 3/2 (— ) 3 ^ 3/2 (^o 3 • 

v eq 

(4.19) 

For a = a_2l0 -2 , (3 = 10, v eq = 0.1 and G/j, = (Gfj,)el0~ 6 , the result becomes 

n G (> Mi, t) ~ 4 • 10- e uz- 3 / 2 (t)v- 3/2 a^(Gi,)l /2 /iJo^Mpc)" 3 , (4.20) 
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which is also plotted in Figure 1 for V{ = 0.25. Similarly, (4.13) is replaced by 

tord(t) ~ n c (M c , t)M 2 6nGt 2 ~ ^va^G^qV^fy- 1 z~ 2 {t) . (4.21) 
For the same parameter choices as above, this becomes 

Sl^t) ~ 10" W 2 (G//)iur 1 *(t)- 2 4 ) . (4.22) 

Eqs. (4.20) and (4.22) are the main results of our calculations. 

As a final consistency check we must verify that M c (t) < M2(t). This is indeed 
true provided that 

z(t) < ^ZeqV^v-^Pa-^Gpi ~ 17/i| (4.23) 

(for Vi = 0.25). However, there may be an even more restrictive condition. Since 
the accretion onto a moving string is not spherically symmetric, M 2 might be 
determined not by (4.3) but by the stronger criterion 

q x (ti,t) > octiz(ti) . (4.24) 

Using (3.26), this becomes 

z(ti) > ^ap-^G^Vizit) = z m . (4.25) 

In order for our results (4.20) and (4.22) to be valid, z m needs to be smaller than 
z(t c ). This will only be the case if 

-PGfJLZegV^V- 1 > z(t) , (4.26) 

which is marginally satisfied for z — 4 if V{ — 0.25 and v eq = 0.1, 

z(t) < 3/i| . (4.27) 

For values of z larger than 4, the values of no{t) and Q n i are suppressed beyond 
the results (4.20) and (4.22) since only the tail of the loop ensemble with velocities 
smaller than the mean velocity V{ = 0.25 manage to accrete a substantial amount 
of mass. 
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5. Discussion 



We have studied the accretion of hot dark matter onto moving cosmic string 
loops and made use of the results to study early structure formation in the cosmic 
string plus HDM model. Our main results are expressed in Eqs. (4.20), (4.22) and 
(4.26). 

The loop accretion mechanism is able to generate nonlinear objects which could 
serve as the hosts of high redshift quasars much earlier than the time cosmic string 
wakes start turning around (which for Gji = 1CT 6 and v s = 1/2 occurs at a 
redshift of about 1). However, there is an upper cutoff to the redshift of large 
mass objects which form by this mechanism given by Eq. (4.26) and for V{ = 0.25 
it corresponds to a redshift of about 4. For larger redshifts, only the loops with 
velocities sufficiently small compared to the mean loop velocity will be able to seed 
nonlinear objects. Note that this redshift cutoff is independent of the parameters 
a and v of the cosmic string scaling distribution which must be obtained from 
numerical simulations and are still quite uncertain. 

The fraction Q n i(z) of the total mass accreted into nonlinear objects by string 
loops unfortunately depends very sensitively on a and v. On the other hand, this 
is not surprising since the power of the loop accretion mechanism depends on the 
number and initial sizes of the loops, and the scaling relation Q n i ~ va is what 
should be expected from physical considerations. 
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Comoving space density of quasars (for h=1/2) 
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Figure 1: Comparison of the observed number density of quasars (+ marks) with 
the number density of protogalaxies predicted in the cosmic string theory with hot 
dark matter for the parameters discussed in the text. The horizontal axis is the 
redshift. 



For the values v — 1 and a_2 = 1 which are indicated by recent cosmic string 
evolution simulations 17 ) , we conclude from Eq. (4.20) that the loop accretion mech- 
anism produces enough large mass protogalaxies to explain the observed abundance 
of z < 4 quasars (see Figure 1). Note that the amplitude of the predicted proto- 
galaxy density curves depends sensitively on the parameters of the cosmic string 
scaling solution which are still poorly determined. Hence, the important result 
is that there are parameters for which the theory predicts a sufficient number of 
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protogalaxies. Since not all protogalaxies will actually host quasars, and since 
the string parameters are still uncertain, it would be wrong to demand that the 
amplitude of the uq curve agree with that of the observed uq. 

It is more difficult to make definite conclusions regarding the abundance of 
damped Lyman alpha absorption systems. In the form of Eq. (1.2), the condition 
for the cosmic string loop accretion mechanism to be able to explain the data is 
also satisfied. However, Eq. (1.2) refers to the value of Q in baryonic matter. The 
corresponding constraint on the total matter callapsed in structures associated 
with damped Lyman alpha systems is 

n BL (z<3)> f^lO' 3 (5.1) 

where fa is the local fraction of the mass in baryons. From Eq. (4.22) it follows that 
the above constraint is only marginally satisfied, and this only if the local baryon 
fraction fa exceeds the average value for the whole Universe of about fa = 0.1. This 
could be another manifestation of the "baryon crisis" for galaxy clusters 28 ), the 
fact that fa in clusters seems to exceed what is expected based on nucleosynthesis 
constraints in a f2 = 1 Universe. On the other hand, in the cosmic string model 
with HDM we expect fa in nonlinear objects to be enhanced over the average fa 
because baryons are able to cluster during the time that the HDM is prevented from 
accreting by the free streaming 29 -*. Thus, cosmic strings may be able to explain 
the baryon excess in clusters and restore agreement with (1.2) in a natural way. 
More calculations are required to resolve this issue. 

In conclusion, we have established that in addition to being in agreement with 
large-scale structure and CMB data, the cosmic string and hot dark matter model 
is also able to produce a sufficient number of protogalaxies at redshifts z < 4 which 
could explain the observed abundance of quasars. A prediction of the model is that 
the distribution of these quasars should be less correlated with today's large-scale 
structure than in inflation-based models, since the loops giving rise to quasars 
are not correlated with the long strings present at the time of equal matter and 
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radiation, which give the dominant contribution to today's large-scale structure. 
However, some correlation might still be present since the loops are correlated with 
the long strings from which they were split off, and since the quasar host galaxies 
evolve into present-day galaxies and may fall into the potential wells created by 
the string wakes. 
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